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Abstract 

In [arXiv:0804.2630], we have analyzed the leading logarithms of energy that appear in 
the inclusive spectrum of gluons produced in heavy ion collisions, calculated in the Color 
Glass Condensate framework. The main result of this paper was that these logarithms are 
intrinsic properties of the colliding projectiles, and that they can be resummed by letting the 
distributions of color sources in the nuclei evolve according to the JIMWLK equation. 

An essential step in the proof of this factorization result is the calculation of the gluon 
spectrum at Next-to-Leading order, and in particular a functional relationship that expresses 
the NLO correction as the action of a certain operator on the LO spectrum. 

In this paper, we show that this type of relation between spectra at LO and NLO is not 
specific to the production of gluons, but that it is in fact generic for inclusive spectra in 
heavy ion collisions. To illustrate this, we compute up to NLO the inclusive spectrum of some 
hypothetical scalar fields, either color neutral or colored, that couple to gluons. 

1 Introduction 

Heavy ion collisions at high energy involve in a crucial way the partons with a very small longitu- 
dinal momentum fraction x of the the nucleon to which they belong. For instance, in heavy ion 
collisions at the LHC, where the center of mass energy is 2.76 TeV per nucleon pair and where 
the typical produced particle has a transverse momentum of only a few GeVs, one probes partons 
with a momentum fraction x ~ 10 -3 or lower. In this region, the valence quarks are completely 
negligible, and most of the partons are in fact gluons (see [1] for instance). The distribution of sea 
quarks-produced by the splittings g — > qq-is suppressed by a power of a s compared to the gluon 
distribution. 

As long as the gluon occupation number remains of order one or smaller, the evolution with x 
of the gluon distributions is governed by the (linear) BFKL equation [2, 3]. The evolution of these 
distributions with x is driven by soft gluon radiation : when one computes inclusive observables 
beyond leading order, they contain contributions that have large logarithms of 1/x. It turns out 
that these logarithms have a universal structure, and can be resummed simply by making the 



gluon distributions ^-dependent (with an x-dependence dictated by the BFKL equation) . This is a 
particular example of factorization, known as £; t -factorization in this case since it involves transverse 
momentum dependent gluon distributions 1 [10-12]. It should be emphasized that causality plays 
a crucial role in factorization: it is because the two projectiles cannot be in causal contact before 
they collide that their partonic content is an intrinsic property that does not depend on what they 
collide with, nor on what observable is going to be measured in the final state. 

In the regime of very small values of x, probed in high energy heavy ion collisions, one may 
reach values of the gluon occupation number that are much larger than one. This is in fact a 
consequence of the BFKL equation, that predicts a growth of the form x~ x , where A is some positive 
exponent. When the occupation number is large, non-linear interactions among the gluons, such as 
recombinations, become important [13]. These processes tend to stop the growth of the occupation 
number when it reaches a value of the order / ~ l/a s , a phenomenon known as gluon saturation. 
The domain where saturation effects are important is delineated by a characteristic momentum 
scale, the saturation scale Q s , that depends both on x and on the nucleus size A (it increases for 
smaller x and for larger nuclei). 

Computing observables in the saturation regime is a challenge for several reasons. Firstly, since 
gluon recombinations are important, processes may be initiated by more than a single gluon from 
each projectile. The consequence is that in order to predict anything, one needs informations 
about the distribution of multigluon states in the wavefunctions of the two incoming projectiles. 
Given these multigluon distributions, the second difficulty is to compute observables from them. 
Indeed, when the gluon occupation number is of order a" 1 , computing an observable at a given 
order in a s requires to resum an infinite set of Feynman graphs. Finally, as was already the case 
in the dilute regime, these observables receive large logarithmic corrections at higher orders. But 
now, these logarithms may be more complicated because their coefficients may be sensitive to the 
multigluon states mentioned above. It would be highly desirable to show that these logarithms arc 
still universal (a property that we expect to be true from causality), and that they can be factorized 
into evolved multigluon distributions. 

The Color Glass Condensate (CGC) is an effective theory, based on QCD, in which these com- 
putations can be carried out (see [14-20] for an account of the early developments and [21-23] for 
reviews). The main simplification compared to full fledged QCD stems from the realization that 
the fast partons in a nucleon or nucleus have a large Lorentz boost factor that slows down their 
dynamics, so that they can be considered as static for the duration of a collision at high energy [24- 
26]. This allows to describe them as time independent color currents that fly along the trajectories 
of the two projectiles (in the light-cone directions z — ±t), whose transverse distribution reflects 
the position in the transverse plane of these fast partons at the time of the collision. All we know 
about these currents is the probability distribution W[p] for the density p a of the color charges that 
produce them. In contrast, no simplification is made to the dynamics of the slower partons, and 
therefore they are described as standard gauge fields, eikonally coupled to the color current of the 
fast partons. The two types of degrees of freedom (color sources and gauge fields) are separated 
by a cutoff A on the longitudinal momentum 2 . In order to compute the expectation value of an 

This form of factorization has been argued to remain valid for some observables (namely, the single inclusive 
gluon spectrum) even when one of the projectiles has a large gluon occupation number, e.g. in proton-nucleus 
collisions [4-9]. But the transverse momentum dependent gluon distribution that describes the nucleus in this case 
is not the expectation value of the number operator. Instead, it has been shown that this distribution is the Fourier 
transform of a two-point correlator of Wilson lines, the so-called dipole amplitude. 

2 In the case of a collision, there is one such cutoff for each nucleus: a cutoff A+ on the momentum fc+ for the 
right moving nucleus and a cutoff A~ for the momentum k~ in the left moving nucleus. 
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observable in this framework, one first computes it for a fixed p, and then averages the result over 
all possible p's according to the probability distribution W[p\. 

In the CGC framework, the power counting appropriate to the description of the saturated 
regime should assume that p ~ g^ 1 , in order to select all the relevant graphs when the gluon 
occupation number is of order l/a s . Under this assumption, one sees easily that observables can 
be expanded in powers of g 2 , each successive order corresponding to an additional loop. The main 
difference with the g 2 expansion in the dilute regime is that there are now an infinite set of graphs 
contributing at each order - the leading order (LO) corresponds to all the tree diagrams, the next- 
to-leading order (NLO) to all the I-loop graphs, etc. Moreover, the tree diagrams that contribute to 
inclusive observables can be organized in such a way that they contain only retarded propagators 
[27]. Consequently, the gluon inclusive spectrum at LO can be expressed in terms of retarded 
classical solutions of the Yang-Mills equations [27]. Likewise, the quark spectrum is given at LO 
by solutions of the Dirac equation in the presence of the previous classical fields [28, 29]. 

Because the CGC is an effective theory with two types of degrees of freedom, loop integrals 
should be cutoff at A^ 1 to avoid double counting contributions that are already included via the 
color currents that describe the fast partons. However, some of these loop integrals contain log- 
arithms of the cutoff - making observables depend on an unphysical parameter that has been 
introduced by hand. In the case of reactions involving a single nucleus, such as Deep Inelastic Scat- 
tering, it is well known that these logarithms can be absorbed into a redefinition of the probability 
distribution This redefinition amounts to letting W become cutoff dependent, with a cutoff 

dependence controlled by the JIMWLK equation. When W[p] obeys the JIMWLK equation, the 
cutoff dependence in W precisely cancels the cutoff dependence from the loop integrals, leading to 
observables that are cutoff independent. Moreover, it is this scale dependence of the distributions 
W that gives observables their rapidity dependence (for instance, the gluon spectrum computed in 
the CGC framework at LO is totally independent of rapidity, and becomes rapidity dependent after 
these logarithms have been resummed). 

Proving the factorization of the logarithms of the cutoffs in the case of nucleus-nucleus collisions 
is more complicated. The main difficulty is that-already at LO-observables cannot be calculated 
analytically 3 because this involves solving the classical Yang-Mills equations in the presence of 
the two strong currents that describe the fast partons of the two nuclei. At NLO, one needs to 
compute I-loop graphs in this classical background field, which again cannot be done analytically. 
This is where causality becomes very handy : the above mentioned technical difficulties arise only 
in the forward light-cone, i.e. after the collision has happened, while the physics that controls 
the logarithms we want to compute happens before the collision, in the wavefunctions of the two 
projectiles. In [37], we have developed a method for separating these two stages of the time evolution 
in the case of the inclusive gluon spectrum, and we have shown that this separation indeed allows 
one to extract analytically the leading logarithms and to prove that they can be factorized in 
JIMWLK evolved W distributions (one for each nucleus). At the moment, this factorization has 
been shown for inclusive observables involving only gluons. For instance, the computation of the 
inclusive quark spectrum done in [28, 29] assumes that such a factorization also works for quark 
production, although this has not been demonstrated so far. 

The cornerstone of this method is a functional relationship between the spectrum at LO and 
at NLO, that makes this separation explicit (see eq. (6) for this formula in the case of the single 
gluon spectrum). It was later found that the same separation can be made for more complicated 
observables, such as the inclusive multigluon spectra [38, 39], or the energy- momentum tensor, 

3 Of course, such computations are routinely doable numerically by now [30-36]. 
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thereby allowing one to prove the factorization of their logarithms of A* . An important question 
toward extending factorization to other observables is the range of validity of these relationships 
between their tree-level and 1-loop contributions. In this paper, we show that such formulae also 
exist in theories that have not only gluons, but also other fields that couple to gluons. For sim- 
plicity, we consider only scalar fields, but the extension of our second example to quarks would be 
straightforward (apart from the complications related to dealing with spinors rather than scalars). 

In the section 2, we provide a brief reminder of existing results concerning the gluon spectrum at 
LO and NLO, and the factorization of its logarithms. More importantly, we provide a diagrammatic 
method for manipulating the objects and operators that arise in the formulae we wish to generalize. 
All these pictorial representations are in a one-to-one correspondence with some equations, but 
provide a much simpler and intuitive way of manipulating them. In the section 3, we extend this 
result to some hypothetical color neutral scalar particle that couples to a pair of gluons. Because 
such a field does not see directly the background field, the results of pure Yang-Mills theory also 
apply here without any change. The main part of the paper is the section 4, where we consider 
colored scalar fields that live in the adjoint representation of the gauge group. These fields have 
nontrivial interactions with the background gauge field. The main result of this section is the 
eq. (60), a natural generalization of the formula we had obtained for the gluon spectrum, that 
relates the LO and NLO spectra. The section 5 is devoted to a summary and concluding remarks. 



2 Reminder on gluon production 
2.1 Gluon spectrum at LO and NLO 

The single inclusive gluon spectrum can be expressed as follows in terms of the color field operator 
A* : 



where p = \p\ is the on-shell energy of the produced gluon and e]^ (p) is the polarization vector. 
In the inclusive spectrum, one should sum over all the possible polarizations A and colors a of the 
produced gluon. The right hand side of this formula involves the two-point H — Green's function 
of the Schwinger-Keldysh formalism 4 . This is reminded by the subscripts + and — carried by the 
two field operators. 

In the Color Glass Condensate framework, gluons are coupled to an external classical color source 
J M that describes the fast partons contained in the two projectiles. This means that there is a non- 
zero disconnected part in the above two-point function, that must be taken into account. Moreover, 
in the regime where the two projectiles are saturated (as is the case in heavy ion collisions), one 
should assume in the power counting that this source is parametrically of order g~ x . Under these 
conditions, the two-point correlator can be expanded as follows 

(Al a (x)A^(y)) = A» a (x)A va (y)+A» a {x)a va {y) + a^(x)A» a (y) + G^ a (x, y) + ■ ■ ■ (2) 

v v ' v- ' 

LO NLO 

The leading order contribution comes entirely from the disconnected part of the two-point function, 
where A^ is the solution of the classical Yang-Mills equations, with null boundary conditions in the 



4 As opposed to a time-ordered Green's function in scattering amplitudes. This peculiarity can be seen for instance 
from the fact that the gluon spectrum is the expectation value of the number operator a^a. 
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remote past, 



[V^F V ] = J v , F v = d"A u - d t/ A" + iglA", A v ] , lim A" = . 



(3) 



From this classical color field, the LO gluon spectrum is given by 



d 3 p 



(2 



^— ^eW(p)eW*(p) J d A xd A y e^"^ D x D y A»(x)A» a {y) 

A, a 



(4) 



The classical field A^ that enters in this formula is not known analytically in the case of a collision 
between two saturated nuclei, but it is fairly straightforward to obtain it numerically [30-36] . Note 
also that, at LO, the inclusive gluon spectrum is of order 



d 3 p 



0(g- 2 ) 



(5) 



when the two projectiles are saturated. 

At next to leading order, the calculation of the gluon spectrum involves small perturbations to 
the above classical field. As was shown in [37], these corrections can be combined into a compact 
functional relationship that expresses the NLO in terms of the LO: 



dN n 



d 3 p 



Jd 3 u [a-T u ]+^Y^J (2 ^ 2fc Jd 3 ud 3 v [a- k xa • T„] [a +kXa ■ T v ] 



dNg 
d 3 p 



(6) 



In this formula, the inclusive spectrum at leading order should be viewed as a functional of the 
classical color field A^ on some Cauchy surface 5 S. The fields a and a±k\ a are small perturbations 
to the classical color field on £, that are defined in detail in [37]. The operator T u is the generator 
of the shifts 6 of the classical field A at the point ueS. 



2.2 Leading logarithms and factorization 

The operator in the square brackets in the right hand side of eq. (6) contains momentum integrals - 
the explicit integration over d 3 k in the second term, and a loop integral hidden in the definition of 
a. These integrals are both divergent when the longitudinal components of the momentum, k + or 
fc~, go to infinity. In the CGC framework, these integrations should have upper limits A ± , where 
A* is the cutoff that separates the fields from the classical sources, in order to prevent double 
countings. The k integrals are thus finite, but they give logarithms of the unphysical cutoffs A*. 

5 This surface S could be any locally space-like surface, such that the field At 1 above S is completely determined 
by the knowledge of its value and that of its first time derivative on S. u is the 3- volume measure on this surface. 

6 The retarded classical color field Af* implicitly depends on its initial condition A s on the initial Cauchy surface 
E. Any functional F[A] can thus be viewed as a functional G[_4 E ] of the initial condition. The value of this functional 
for the new classical field resulting from a shift of the initial condition can be obtained by exponentiating the operator 
T„, 

G[A S + a] = cxp [y d A u [a ■ T„]] G[A S ] . 

E 

(this relationship can be seen as a definition of T„.) 
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In [37], the following result was established 

d 3 k 



J d 3 u [a-T u ] + l -Y,J ( 2 ^)32fc / [ a - kXa ' [° +feAa ' T "] = 

= log(A+) Hi + log(A") n 2 + terms w/o logs , (7) 

where the operators Hi and 7^2 are the JIMWLK Hamiltonians of the two nuclei. The most 
important aspect of this formula is that the logarithms of the two cutoffs are multiplied by objects 
that depend only on the sources of the corresponding nucleus. This property is crucial for the 
universality -and thus the factorization- of these logarithms: it ensures that the logarithms of A + 
are an intrinsic property of the projectile moving in the +z direction, that does not depend in any 
way on the projectile moving in the opposite direction. 

Thanks to eq. (7), one can resum all the leading powers of the logarithms (i.e. all the terms of 
the form (a s log A 1 * 1 )™), and absorb them into the scale dependence of the distributions Wi[pi] and 
W2[p2\ that describe the source content of the two projectiles: 



dNg 
d 3 p 



= J [D Pl D P2 ] W x \px]W 2 [p2] 



dN g 



Leading Log J a P 



(8) 



where the two distributions W obey the JIMWLK equation 



A ± ^^ = -W 1 , a Wi, 2 . (9) 



2.3 Diagrammatic manipulations with T u 

In the above NLO expression (6), a very compact result is obtained thanks to the introduction 
of the operator T u . This operator will play a crucial role in the rest of this paper, and for this 
reason it is very useful to gain some intuition on how it acts on various objects. All the formulae 
involving T u can be proven by making use of the Green's formulae that relate classical fields - 
and perturbations thereof- to their initial value on the surface E. However, these proofs are often 
cumbersome (see [37, 38] for some examples). In this subsection, we would like to present a more 
intuitive, diagrammatic, way of manipulating these operators. 

First of all, let us recall the most basic identity (derived in [37]), that relates a linearized 
perturbation propagating over the classical background field A* 1 to the classical field itself 

a" (a:) = J d 3 u [a ■ T u ] A"(x) . (10) 

This formula means that if we know the perturbation on the surface S and how the classical field 
A(x) depends on its initial condition on E, then we can obtain the value of the perturbation at 
the point x by acting on A(x) with the operator TT U . Roughly speaking, T u acts as a first order 
derivative with respect to the initial value of A (and that of its first time derivative) at the point 
u E E. Diagrammatically, the identity of eq. (10), before the integration over d 3 u is performed, 
can be represented as follows : 




G 



In the left hand side, the thin wavy line terminated by a gray blob represents the retarded classical 
color field at the point x. In the right hand side, one has a perturbation that starts on E as 
a(u), and then propagates to the point x over the classical background (the boldface propagator 
thus indicates that it is dressed by this background field). In this diagrammatic representation of 
eq. (10), one should keep in mind the following rule : when the combination [a ■ T u ] acts on some 
object, whatever is attached to the operator T u (here a(u) ) replaces the part of that object that 
hangs below S at the point ueS. 

Let us now consider the action of an operator that contains several powers of T, the simplest of 
which is [a • T u ][6 • T„], that involves two T's at different points u, v £ S. When this combination 
acts on the retarded classical field A 11 , one gets an object that involves the dressed 3-point vertex 
with a(u) and b(v) attached to two of its endpoints: 



[a-T u ][b-T v ] 



(12) 



a(u) b'y) 

In the right hand side of this relationship, the blob indicates that the 3-gluon vertex is dressed by 
the background field, 

(13) 

i.e. it is the sum of the bare 3-gluon vertex and a 4-gluon vertex where one of the four legs is 
attached to the classical background field. 

In order to understand diagrammatically the equation (6), we need also the following represen- 
tations, 




E 



d 3 k 
{2n) 3 2k 



a(u) 



a-kXa(u)a+kXa(v) 




(14) 



The first one is just the definition of a(u) as the 1-loop correction to the field expectation value 

at the point u £ S. The second equation is a representation of the G-| propagator (dressed 

by the background field) as a sum over a complete basis of fluctuations that propagate over this 
background (see [37], and also eq. (40) later in this paper). It is now easy to stitch eqs. (11), (12) 
and (14), in order to get 



d 3 U 




d 3 u d 3 v 




(15) 
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The two terms in the left hand side of the second line correspond to the two possible localizations 
of the 3-gluon vertex: below or above 7 the surface E. When combined, these two terms reconstruct 
the complete 1-loop correction to the expectation value of the field at the point x (the prefactor 
1/2 in the second term gives the correct symmetry factor for the tadpole). 

Since T u acts as a derivative, its action on a product of classical fields is distributed according to 
Leibnitz's rule. Since the gluon spectrum at LO is the Fourier transform of the product A(x)A(y), 
it is instructive to see how T u acts on this product. For instance, on gets 




a(u) b(v) a(uj b(v) 

Combining these two equations with eqs. (14), one sees graphically that eq. (6) generates exactly 
the terms that are needed to obtain the gluon spectrum at NLO, 




This diagrammatic approach is very powerful, since the above manipulations are arguably much 
simpler than the derivation of the NLO gluon spectrum in [37] . Alternatively, if one does not go all 
the way to relying only on these diagrams, one can use them as a guidance to guess some identities, 
that can then be proven by the more analytical (and tedious) methods of [37]. 



3 Color neutral scalar particles 
3.1 Model 

The simplest extension to the pure Yang-Mills case is to consider the production of some hypothet- 
ical scalar particle, singlet under the gauge group of QCD, that couples to the square of the gluon 
field strength. Such a particle can be described by a field 4> whose interactions are given in the 

7 A crucial aspect in all these manipulations is the fact that the inclusive quantities can always be expressed 
in terms of retarded propagators, that describe the causal propagation of some object over a background field. 
Therefore, the "stitching" procedure described above simply corresponds to concatenating the evolution from — oo 
to the surface E (encoded in the prefactors of the T's) and the evolution from S to the time x° of interest (encoded 
in the object upon which the T's act). 
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following Lagrangian density, 

r=I(5^)(^)_l m 2^_ ^g^F^p^a (ig) 

v v ' v v ' 

free scalar field coupling to gluons 

For instance, this particle could be a Higgs boson, and this Lagrangian would be an effective 
description of its coupling to a pair of gluons. At the fundamental level, this coupling occurs 
mostly via a top quark loop, and one would obtain this Lagrangian by integrating out the quarks, 
leading to an effective coupling constant g^. Note that this model was considered recently in [40], in 
the limit where the scalar particle mass is much larger than the saturation momentum. This limit 
is simpler because the leading term (when we expand both in g 2 and Q 2 /m 2 ) in the production 
amplitude has one attached to each of the colliding nuclei, with no entanglement between 
the sources of the two nuclei, leading to a much simpler form of factorization. In this section, we 
consider a generic mass m 2 for this scalar particle, possibly as small as the saturation scale, and 
therefore we do not expand observables in Q 2 /m 2 . 

Here, we disregard the self-interactions of these scalar particles, and consider only their coupling 
to gluons. Moreover, in the power counting, we assume that the effective coupling constant g^ is 
much smaller than the strong coupling constant g. Therefore, we only consider contributions that 
have the lowest possible order in g^, and the perturbative expansion is in powers of g. This means 
that the NLO corrections affect the gluons, and only indirectly the scalar field 4> via its coupling to 
F 2 . 

In this model, one can see the gluon field strength squared, multiplied by the coupling g^, as a 
source for the scalar <j>, 

1 

The Lagrangian (19) is thus that of a free scalar field coupled to the external source J. 



J=— Ag<j>F ; v F^ a . (20) 



3.2 Inclusive spectrum at Leading Order 

In close analogy with the gluon spectrum (1), the single inclusive spectrum for cf> is given by the 
following formula, 



d 3 p (2n) 3 2E p 



d 4 xd 4 ye tp - {x -y> (n x + m 2 )(n y + m 2 )(Mx)4>-(y)) , (21) 



in terms of the two-point H — function for the field <fi- This formula is true to all orders, but at 
leading order in g 2 the two-point correlator that appears under the integral is equal to the product 
ip(x)ip(y), where ip(x) is the retarded solution of the classical equation of motion for <fi, 



{U x + m 2 ) V >{x) = --g^ v F^ a , lim V (x) = . (22) 



In this equation, it is sufficient to keep only the leading order contribution to the gluon field 
strength, i.e. its classical value given by eq. (3). In the regime where the colliding projectiles have 
a saturated gluon content, the gluon field strength is of order 1/g. Thus the source J is of order 
g^/g 2 and the leading contribution to the <j> spectrum is of order g^/g . In contrast, the first 
connected contribution to (c/> + (x)<f)-(y)) starts only at the order g^/g 2 - 
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Qs 



V 



x> cf 



Figure 1: Diagrammatic representation of the contributions to the inclusive <f> spectrum at Leading 
Order. The dotted line is the free propagator of the scalar particle (eliminated in the spectrum by 
the action of the operators □ + m 2 ). Each spiraling line terminated by a circular blob represents 
one instance of the gluon field strength at tree level (i.e. it contains terms of order one and two in 
the gauge field A* 1 ). 

Thus, at LO, the </> spectrum is a sum of tree diagrams (see the figure 1). In terms of the 
classical field strength, the inclusive spectrum for (j) at LO reads simply 



dN* 



d 3 p 



1 



(2n) 3 2E p 16 



d 4 xd 4 y e^-s/) jra {x)F ^a {x) {y)F P^ {y) _ (23) 



This formula simply expresses the well known fact that, when non-self- interacting particles are 
produced by an external source, the inclusive spectrum is the modulus squared of the Fourier 
transform of the source. 

3.3 Next-to-Leading Order correction 

Let us now consider the production of the ^-particles at NLO, i.e. at the next order in the strong 
coupling constant g 2 . At the next order, one has two types of contributions: 

i. 1-loop corrections to the gluon field strength, or to the gluon field strength squared. 

ii. the connected part of ((f> + (x)<f)_(y)) at tree level. 

These two contributions are illustrated in the figure 2. These corrections affect only the gluons, and 



X 




y x f™\ °>\^ y x 



Figure 2: Diagrammatic representation of the contributions to the inclusive <j) spectrum at Next- 
to-Leading Order. Each blob represents an infinite sum of tree level graphs. The diagram on the 
left contains the 1-loop correction to F^ v (x). The diagram in the middle involves the connected 
part of the 1-loop correction to F ilu {x)F )lv {x) (there are two other topologies, not represented here, 
where the 1-loop correction is on $_(?/)). The right diagram contains the lowest order connected 
contribution to F^ v (x)F pa (y). Thick lines indicate that the gluon propagators are dressed by the 
classical background field 

they can be handled in exactly the same way as the corrections to the inclusive gluon spectrum. 
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The manipulations performed in the previous section for the expectation value of A^(x)A v {y) 
can be extended to this more complicated situation, leading the following formal expression 8 for 
the inclusive <fi spectrum at NLO, 



dN* 



d 3 p 



fd 3 u [a-T u ] 



L £ 



2 

A, a 



d 3 k 
(27r) 3 2fc 



J d 3 ud 3 v [a^ k xa ■ T u ] [a +kXa ■ T„] 

y. 



dNs 



d 3 p 



(24) 



In this formula, the inclusive spectrum at leading order should again be considered as a functional 
of the classical color field on some Cauchy surface S. This relationship is formally identical to 
the one that relates purely gluonic observables at LO and NLO (e.g. inclusive gluon spectra, or the 
gluon contribution to the energy-momentum tensor). 

In the gluon case, we have shown that it is the integration over k in the right hand side of eq. (24) 
that leads to the logarithms of the CGC cutoffs. Since we have here exactly the same operator, the 
logarithms that appear in the NLO correction to the inclusive <$> spectrum are identical to those 
encountered in the gluon spectrum, and they can be resummed into the JIMWLK evolution of the 
distributions of the sources that produce the color field, 



d 3 p 



Leading Log 



= J [D Pl D P2 ] Wi\pi]W 2 \p2 



d 3 p 



(25) 



This is another example of the universality of these logarithms. 



4 Adjoint scalar s 
4.1 Model 

The simplicity of the previous example is to a large extent due to the fact that the field <j> is not 
charged under the gauge group of strong interactions. Therefore, the NLO corrections are entirely 
due to corrections to the gauge fields themselves, and to their two-point correlations. 

Things get more complicated when the produced particle carries a color charge. An obvious 
example is of course that of quarks and antiquarks. However, since quarks also bring the compli- 
cation of being described by anticommuting fermion fields, let us consider as an intermediate step 
the case of some hypothetical colored scalar field 9 a , that lives in the adjoint representation of the 
gauge group. The Lagrangian for such scalars is 

L^(D,9) a (D^Y a , (26) 

where D„ is the covariant derivative that ensures the minimal coupling of the scalars to the gluons. 
We assume here that these scalar fields do not interact directly among themselves (i.e. their 
interactions are always mediated by gluons). 

8 One can easily see that this derivation applies to any quantity that can be expressed at leading order as a local 
function of the retarded classical color field. It also applies to bi-local expressions —i.e. involving the gauge field at 
two points x and y— provided one of these points carries the Schwinger-Keldysh index + and the other one the index 
— . For more than two points or more general assignments of the Schwinger-Keldysh indices, it still works provided 
all the pairs of points have space-like separations. 



11 



4.2 Inclusive spectrum at LO 



The starting point to compute the single inclusive spectrum for the 9 particles is the analogue of 
eq. (21), 

S| = (2^E / ^V- (I - S » □.□„ (e%{xW_{y)) , (27) 

where c is the color of the produced scalar particle 9 (summed over in the definition of the inclusive 
spectrum). The two-point correlator under the integral differs from the previous two examples in 
the fact that it has only connected contributions. This is due to the fact that the Lagrangian (26) 
is even in the field 8. To stress this property, let us rewrite eq. (27) as follows 

f| = Wf2 P ^ j d^y U x U y D-_(x,y) , (28) 

where D b £_ is a connected two-point function. The Leading Order contribution to D+_, that we 
will denote £>+_, is the sum of all the tree level graphs, i.e. the two-point function dressed by the 
classical color field _4 M , solution of eq. (3). Tr£_ can be expressed in terms of a complete basis of 



Figure 3: Diagrammatic representation of the only contribution to the inclusive 6 spectrum at 
Leading Order (disconnected graphs cannot exist because the Lagrangian is even in 6). The thick 
line indicate that the propagator is dressed by insertions of the classical color field Ab- 
solutions of the scalar wave equation on top of the classical color field A^, 

{V^) cb d b ±qa {x) = 0, lim # b ±qa (x) = date*** . (29) 

In the notation $± qa (x), one should not confuse the two color indices a and b. The lower index 
a indicates the color of the wave in the remote past, before it propagates over the classical color 
field, while the upper index b indicates the color at the point x. The fields "&± qa can be represented 
diagrammatically as follows : 

x 



$± qa (x)= X . (30) 




)q,a 

In this representation, the red cross represents the initial condition at a: — > — oo, i.e. a plane wave 
with fixed momentum and color. The arrow is a reminder of the fact that the functions $±q a obey 



9 In the following, when the context is not ambiguous, we may omit the color indices attached to the endpoints of 
the objects we are manipulating, in order to keep the notations lighter. 
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retarded boundary conditions set at x° — > — oo. In other words, the propagator that connects the 
red cross to the point a; is a retarded propagator (dressed by the classical color field A^). 
In terms of these objects, one can represent the LO propagator X>_| as 




(31) 



and the 8 spectrum at LO is given by 



dNg 

d 3 p 



(2 



(32) 



In this formula for the 9 spectrum at leading order, q can be interpreted as the momentum of 
the antiparticle that must necessarily be produced along with the particle of momentum p that is 
tagged in the measurement. 



4.3 NLO corrections 

4.3.1 List of the 1-loop corrections 




Figure 4: Diagrammatic representation of the contributions to the inclusive 9 spectrum at Next-to- 
Leading Order. All the lines represent propagators dressed by insertions of the classical color field 
A^ . The blobs indicate that the corresponding vertices are dressed by the background field. 

The NLO corrections to the inclusive spectrum of the 9 particles are represented in the figure 10 
4. They can be divided in two distinct categories: 

i. Tadpole insertions on the propagator of the 9 field, that appear in the first two diagrams of 
the figure 4. The leftmost correction is the 1-loop correction a M to the classical color field 
A^ that enters in the right hand side of eq. (6). The second diagram is a local self-energy 
correction on the propagator of the 8 field. 

ii. Non-local 1-loop self-energy correction to the scalar propagator, represented in the third 
diagram of the figure 4. 

10 Note that the gg09* vertex is not dressed at tree level. 
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4.3.2 Alternate representation 

Our goal for the rest of this paper is to write the NLO corrections to the 6 spectrum in a form 
similar to eqs. (6) and (24). This will be done in several steps. In this subsection, we first rewrite 
these corrections in a form that resembles eq. (31), in order to highlight the terms that can be seen 
as mere corrections to one of the functions "d± qa - 

This is very simple in the case of the two terms on the left of the figure 4. Indeed, in these two 
terms the self-energy insertion is a tadpole, that by definition corrects the propagator at a single 
point. This implies that one can see them as a correction to one of the two factors "d± qa - Let us 
denote C±qa this 0{g 2 ) correction to l d± qa , 





C± qa {x) = J\ + A (33) 
)q,a <&q,a 

In terms of C±qa, the corresponding NLO corrections to the propagator _D_| reads 

D { :l(x,y) =J2j [#- qa (xX +qa (y) + (- qa (x)# +qa (y)] . (34) 

This equation simply states that each of the two factors "d± qa in eq. (29) has to be corrected in 
turn. Diagrammatically, this can be represented as 

x y x y 

D l °l(x,y)= \ / + \ / + [x<*y], r.K) 





where the unwritten terms, hidden under the abbreviation [x 4-> y], are the two terms where the 
tadpole is attached to the leg that ends at the point y. 

Let us now consider the third diagram on the right of the figure 4. This term is the H — 
scalar propagator, corrected by a non-local 1-loop self-energy insertion. Because the self-energy is 
a genuine two-point function, this correction cannot be seen as a correction of either i9+q Q or $_q a 
in eq. (29). Instead, we must go back to the Schwinger-Keldysh formalism in order to evaluate this 
term. 

Let us call £ the sclf-cncrgy. Summing over all the possible assignments of the internal ± indices 
in the Schwinger-Keldysh formalism, we can first write 

D^l = £>++£++£>+_ - - £>+_£_+£>+_ + £>+_£__£>__ , (36) 

where the T>±± are the four components of the scalar propagator at tree level. To keep the notations 
simple, we have not written explicitly the space-time integrations involved in concatenating the 
self-energy and the tree-level propagators in this equation, but we have of course included the 
appropriate signs to keep track of the fact that the vertices of type — are opposite to those of 
type +. After some simple algebra, one can rearrange the four terms in the right hand side of the 
previous equation as follows 

£><?>. = 2> R £ fl P+- + T> B V + JD A + V + _V A V A , (37) 
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where the retarded and advanced propagators and self-energies are related to the Schwinger-Keldysh 
ones by 





= v ++ 








v A 


= v ++ 




= v+- 




£„. 


= s ++ 


-£+- 


= £-+ 






= s ++ 


-£-+ 


= x+- 





(38) 

It is interesting to note the pattern 11 in eq. (37) : each term in this formula contains exactly one 
object of type H — . Everything on the left of this object is retarded, and everything on its right is 
advanced. The three terms of eq. (37) exhaust all the possible locations of the object that carries 
the H — indices. 

For the factors X>_| , we can use eq. (29). In coordinate space, the 1-loop self-energy E_| is the 

product of a scalar and a gluon propagator : 



(39) 



where the symbol oc means that we have not written the vertices (they are a bit cumbersome to 
write because of the possibility of attaching the background field A? to one leg of the ggdO* vertex 
to form an effective gOO* vertex - how the complete vertices arise will become clear later). In this 

formula is the gluon propagator dressed by insertions of the classical gauge field . There is 

for a formula identical to eq. (29), 



G+-{x,y) 



E 

a . A 



d 3 k 
(27r) 3 2fc 



a-k\a(x)a+k\a(y) 



(40) 



where the a±k\ a are now gluonic fluctuations that propagate on top of A^, starting as plane waves 
of momentum fc, polarization A and color a when x° — > — oo (see [37] for more details). Mimicking 
eq. (31), we represent this diagrammatically as follows : 



G+-{x,y) 




(41) 



The three terms of eq. (37) can therefore be represented in the following way, 

x y x y 






(42) 



1 This pattern would remain true even when more than three propagators and self-energies are concatenated. 
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where the big grey blobs represent the self-energies T, R A . 

At this point, it is convenient to rewrite the retarded self-energy 12 First of all, one can 
write 

Y, R = E ++ - £ + _ 

oc Q++V++ - Q+-V+- 

= g R v ++ + g + ^v R 

= g R v_ + + g+_v R . (43) 

In order to obtain the last equality, we must recall that 

V ++ (x, y) = d(x°-y°)V„ + (x,y) + 0(y° - x°) V + _ (x, y) . (44) 

Since this propagator is multiplied by the retarded propagator g R , the terms with 8(y° — x°) is 
forbidden, and we can in fact replace the £>++ by T> [_ to get the final formula. In this rearrange- 
ment, we have chosen to break the symmetry between the T> and g propagators in a certain way. 
It is possible to obtain a more symmetric formula. By doing slightly different transformations, we 
could alternatively obtain 

z R <xv R g_ + + v + _g R , (45) 

and, by combining eqs. (43) and (45), 

E a «ix> a (e-+ + 5+-) + i(X> + _+2?_+)e jl . (46) 

Diagrammatically, the expression given in eq. (46) for the retarded self-energy can be represented 
as follows, 




This representation also highlights the different physical meanings of the two terms: 

• In the first term, a vacuum fluctuation made of a pair of virtual gluons is created in the past 
(at the location of the green cross). One of these gluons is absorbed by the scalar fluctuation 
at the point u, and the second of these gluons is absorbed by the same scalar a bit later, at 
the point v. 

• In the second term, the vacuum fluctuation created at the location of the red cross is a 98* 
fluctuation. One of the scalars from this vacuum fluctuation annihilates the incoming scalar 
at the point u to form a gluon. Later, at the point v, this gluon is absorbed by the second 
scalar from the vacuum fluctuation to give the outgoing line. 

12 A similar reorganization of the terms can be performed on the advanced one. 
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4.4 Expression of the NLO corrections in terms of the operators T 
4.4.1 Introduction 

Our goal is now to relate the NLO corrections (at the exception of those that involve 99* vacuum 
fluctuations - see later) to the LO spectrum by means of a formula similar to eqs. (6) and (24). In 
order to keep the derivation as intuitive as possible, we will perform this graphically, by using the 
rules discussed in the section 2.3. Before we proceed, let us summarize the previous discussion by 
listing all the terms involved in the 9 spectrum at NLO : 




[a] [b] [c] [d] 




[e] ffl [g] [h] [i] 



In these graphs, the red and green crosses are objects that are pinned at x° — — oo (because they 
represent the initial condition in the remote past for the scalar waves 'd± qa and gluonic waves a±k\ a ). 
The other vertices in these diagrams can be located anywhere in space-time, and the corresponding 
formulae contain integrals over the positions of all these vertices. In particular, these vertices can 
be located above or below the Cauchy surface £. The only constraint on their locations is that 
they are ordered by the causal property of the retarded propagators: the cndpoint of a retarded 
propagator must lie in the future light-cone of its starting point. 

4.4.2 Action of the gluonic shifts 

As a first step, let us determine the action on the LO spectrum of the operator 




that appears in eqs. (6) and (24). The T's contained in this operator act only on the initial condition 
on E for the gauge field (implicitly hidden in all the propagators, since they are all dressed by the 
classical background field). Let us start with the action of [a-T u ] on the LO spectrum. Graphically, 
we have 



x y x y 



[a-T u ] V + _{x,y) = 
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It is equally straightforward to write the action of the quadratic operator [a ■ T u ][b ■ T„] on the 
leading order expression, 



[a-T u ][6-T„] V+-(x,y) = 




afu) bly) hjv) a) 



When we contract a(u) into T u and § ^2 X a J, a^k\ a (u)a + k\a(v) into T„T„, and integrate over 
u,v <E E, we generate a good deal of the NLO terms in eq. (48). However, by inspection of the 
graphs that are generated in this operation and by comparing with the list in eq. (48), we readily 
see that two types of terms are missing: 

i. The graphs [f,h], where the NLO correction is due to a 89* vacuum fluctuation, are missing 
entirely. Having in mind to use this NLO calculation in order to extract the leading logs of 
energy, we are going to disregard these terms from now on. Indeed, as in the derivation of 
the BFKL equation, only the gluonic vacuum fluctuations contribute at leading logarithmic 
accuracy. Scalar or fermionic fluctuations do not contribute to the logarithms at this order, 
and become necessary only when one goes at the next to leading logarithmic accuracy. 

ii. In the graphs [a,b,c,d,e,g,i], only the portion of the integration domain where the vertices 
u and v are both above the surface E are generated. The missing pieces of these graphs are 
important even at leading log, and in the rest of this section we will show how to generate 
them. 

4.4.3 Scalar shifts 

In the terms of eq. (48), when one of the vertices u or v (or both) is located below the surface S, it 
means that the scalar waves r d± qa have been affected by the gluonic fluctuation before they reach 
E. This implies that the corresponding modification to the final observable should be generated by 
operators that act on the value of the scalar fields on E, instead of the operator T u that acts on 
the gluon fields on E. 

In order to see this, let us introduce the analogue of the operator T u for the fields ■Q± qa - Consider 
a generic solution of the following equation of motion 

(V^) cb #\x) = , (52) 
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that describes the propagator of a scalar over some gauge background (The gauge field in the 
background being fixed.) The Green's formula for the solutions of this equation reads 



#°(x) = i / dyD R (x,y) 2igA, 1 (y)d^ + igid^A^y)) + g A tJl (y)A^(y) ® c {y) 



be 



+i / d 3 uD° R (x,u) (n- du)$ b (u) , 



(53) 



where D R is the free retarded propagator for the scalar field 9 (i.e. the retarded Green's function 
of the D'Alembertian operator D x ). The domain f2 for the integration in the first line is the region 
of space-time located above the surface S. The term in the second line is a boundary term that 
contains all the dependence of the solution on its initial condition on E (n M is the normalized vector 
normal to S at the point u). 

Let us now define the following operator, 



5d b (u) 



(n ■ dt] b {u)) 



8{n-d& b {u)) ' 



(54) 



that mimics for scalar fields the operator T u that was introduced for gauge fields. Diagrammatically, 
the action of this operator can be pictured as follows, 




x 

\ 



(55) 



In words, £ ■ replaces at the point u E Y, the initial value of $ by £(u). Thanks to this operator, 
the perturbation to $ b (x) that results from a perturbation £ b (it) of its initial condition on the 
Cauchy surface E can be written as 



(56) 



When applying this formalism to the spectrum of the 9 field, we must remember that this 
spectrum does not depend on a single $ function, but on a an infinite collection of such functions, 
the $±q a , for all q and a. Therefore, we need an operator such as the one defined in eq. (54) for 
each of them. To avoid the proliferation of notations, we will continue to use the compact notation 
£ • to denote the sum of these operators for all the relevant fields, 



a e=i ' 



(57) 



In order to generate the missing terms in eq. (48) (at the exception of those corresponding to 
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scalar vacuum fluctuations), let us introduce the following two objects 

u u u 



p{u) 



CfcAa(u) 




(58) 



It is now possible to generate the missing terms as follows : 

i. The operator f3 ■ will generate the terms of eq. (48) where the two vertices u and v are 
attached to the same scalar line and are both below the surface S. 

ii. The operator j^C-fcAa ■ [C+^a • I" 5 ] generates the part of the term [i] of eq. (48), where 
the vertices u and v are both below E. 

iii. In order to produce the part of the terms [e,g,i] where one of the vertices u and v is above £ 
while the other is below, we must act with the operator 



1 



1 



a-kXa ■ T][( +kXa • T # ] + -[C-kXa ■ T^a+kXa ■ T] 



(59) 



that mixes the gluonic and the scalar shifts. 
4.4.4 Final form of the NLO spectrum 

Let us now summarize this section by collecting all the terms together. At the exception of the 
contribution of the scalar vacuum fluctuations, the NLO correction to the 9 inclusive spectrum can 
be written as follows, 



dNg 
d 3 p 



Jd 3 u [a ■ T u + /? ■ T£] + 



■E 

i r (ftk r 

+ 2^2J ( 27r )32fc J [ a ~ feAa ' Tu + C - kXa ' T ^ [ Q+fcAa ' Tv + C+kXa ' 



dNg 

d 3 p 



(60) 



This formula is the main result of this paper. It generalizes the eqs. (6) and (24) to the case of a 
theory involving new colored fields in addition to the gluons. The general structure of the operator 
that relates the NLO correction to the LO is preserved by this new formula, and the generalization 
simply amounts to introducing new terms that involve derivatives with respect to the initial value 
of the scalar fields. 

Although we have derived it here for scalar fields, a structurally identical formula can be derived 
for quark production. The NLO topologies in eq. (48) are the same for the quark spectrum, at the 
exception of [b, d] that do not exist for quarks due to the absence of qqgg coupling. We have 
sketched the derivation of the quark analogue of eq. (60) in the appendix A. 



20 



4.5 Discussion on factorization 



In the case of the gluon spectra, the identity (6) was essential in order to extract the logarithms 
of the CGC cutoff A, and to prove that these logarithms can be factorized into distributions of 
color sources that evolve according to the JIMWLK equation. The reason for this is that, if one 
chooses the Cauchy surface E to be just above the past light-cone, then all these logarithms are 
contained in the operator that relates LO and NLO, as stated by the equation (7). From this 
formula, factorization follows easily thanks to the Hermiticity of the JIMWLK Hamiltonian (see 
[37] for details). 

Obviously, the operator in the left hand side of eq. (7) is contained in the operator that appears 
in the right hand side of eq. (60), and as a consequence we already know that the NLO 9 spectrum 
contains logarithms proportional to the JIMWLK Hamiltonian. However, it is crucial to note that 
we do not get the complete Hamiltonian here: we get the JIMWLK Hamiltonian restricted to act 
only on the color sources contained in the gluon fields on S. This is because it is derived from the 
operators T, that act only on the initial gauge fields, but not on the scalar fields. In other words, 
it would be more correct to write 




= log(A+) + log(A") H$ + terms w/o logs , (61) 

where Hf 2 is this restricted JIMWLK Hamiltonian that acts only on gluons. 

For factorization to be valid for scalar production, we would also need the missing parts of the 
Hamiltonian, that act on the initial scalar fields. This is presumably the role of the additional 
terms in the operator that appears in eq. (60). Indeed, these terms involve the new shift operator 
TT^ , that precisely acts on the initial value of the scalar fields. We are thus led to conjecture an 
extension of cq. (7), that reads 

J d 3 u [a ■ T u + (3 • T*] + 

s 

i r d?k r 

+ 2^7 (2 ^3 2fc J [°-*Aa • T u + C-fcAa • T*] [a +kXa ■ T v + C+fcAa • T*] = 

= log(A+)Ki+log(A-)K 2 + terms w/o logs , (62) 

where now Hi. 2 is the full JIMWLK Hamiltonian that acts both on the color sources contained in 
gauge fields and in scalar fields. This formula, if true, would immediately imply the factorizability 
of the leading logarithms of the 9 spectrum into the universal distributions W. However, proving 
eq. (62) requires extracting the logarithms coming from the new terms of eq. (60) and showing that 
their coefficients recombine nicely to form the kernel of the JIMWLK Hamiltonian. We stress that 
this has not been done at the moment. 

5 Summary and outlook 

In this paper, we have derived generalizations of an identity that plays a crucial role in the proof 
of factorization of high energy logarithms in heavy ion collisions. This functional relationship, that 
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relates the leading order and next-to-leading order contributions to inclusive observables, is a central 
step in factorization because it highlights the causal aspects of a collision process, by separating 
what happens before the collision from what happens after. 

The generalization we have obtained goes in the direction of extending the field content of 
the theory under consideration. To the gluons already present in all the known examples of this 
factorization, we have added scalar fields that couple in some way to the gluons. The simplest 
extension is that of color neutral scalar particles. Because they are not colored, the formula that 
was already known for gluons is valid without change in this situation. 

Then, we have considered a much more complicated example, involving colored scalar particles, 
that interact non trivially with the gauge fields. In this case, the formula that was valid for the 
pure glue case must be slightly modified. We have shown that this can be done by a very natural 
extension of the formalism: we also need operators that differentiate with respect to the initial 
value of the scalar fields. Thanks to this generalization, it is also possible to obtain a functional 
relationship between the LO and NLO contributions. 

An obvious extension of this work is to use eq. (60) in order to extract the logarithms of the 
CGC cutoffs, and to prove that they can be absorbed into the JIMWLK evolved W distributions 
that already appear in purely gluonic observables, providing further evidence of their universality. 
Naturally, this would be most interesting in the case of the quark inclusive spectrum, for which we 
have outlined the relevant NLO formula in the appendix A. 

A more remote application of the same ideas and tools is to study the logarithms that arise in 
gluon production beyond leading logarithmic accuracy, in order to obtain the first a s correction 
to the JIMWLK Hamiltonian 13 . The first step in this program is to derive the inclusive gluon 
spectrum at NNLO, i.e. at two loops. Some of the topologies one would need to evaluate in this 
calculation are very similar to the graphs we have studied in the section 4. If one wants to extend 
the factorization proof of [37] beyond leading logarithmic accuracy, an important intermediate step 
would be to find a functional relationship between the gluon spectra at NNLO and LO, similar to 
eq. (6). 

There is another, much more intriguing and difficult, direction in which one could try to extend 
these factorization results - that of exclusive observables, e.g. the cross section for events with a 
rapidity gap. In the case of exclusive observables, the strategy for establishing factorization that we 
have employed in [37] fails at the first stage: so far we have not been able to obtain formulae such as 
eqs. (6), (24) and (60) for any exclusive observable. The main obstacle in obtaining similar formulae 
in this case seems to be that these exclusive observables are expressible in terms of classical fields 
with non-retarded boundary conditions [45] (for instance boundary conditions that constrain the 
fields both at x° = — oo and at x° — +oo). 
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A Inclusive quark spectrum 



The calculation of the inclusive spectrum at LO and NLO for charged scalar fields, that we have 
presented in the section 4, can easily be generalized to the case of quark production. Now the 
Lagrangian density that couples the quarks and antiquarks to the gluons is 

C = V , (63) 

where ip is the spinor for one quark flavor (we disregard the quark mass, that does not play any 
role in this discussion). This Lagrangian is quadratic in the fermion fields, which means that the 
quarks do not interact directly among themselves. In terms of its interactions, this theory is in fact 
a bit simpler than the theory of charged scalars, since it has only a qqg vertex and no derivative 
couplings. 

The inclusive quark spectrum is given to all orders by a formula which is very similar to eq. (28), 
= j^- p E / d ' xdi V eiP<X ~ V) (us(p)0 x )(ui(p)0 y ) S*L(s, y) , (64) 

where u s (p) is a free spinor representing a quark of momentum p and spin s, and S l +_(x, y) is the 
H — component of the quark two-point function in the Schwinger-Keldysh formalism (we use latin 
letters k for colors indices in the fundamental representation). Note that the inclusive spectrum 
is summed over all spin states s and colors i of the produced quark. 

At leading order, we can represent this two-point function in terms of a complete basis of 
solutions of the Dirac equation in the background field A**, 



s . k 

= °> lim Ksk( x )= 6 ^s(q)e iq - X . (65) 

Here also, the lower color index is the color of the quark at x° = — oo, while i is its color at the 
current time a; . Eqs. (64) and the representation (65) were the basis of the numerical evaluation 
of the quark spectrum in the CGC framework done in [28, 29]. 

At NLO, we need to evaluate the 1-loop corrections to the two-point function S-\ The relevant 

topologies are those listed in the figure 5. One of the topologies present in the scalar case does 



Figure 5: Diagrammatic representation of the contributions to the inclusive quark spectrum at 
Next-to-Leading Order. All the lines represent propagators dressed by insertions of the classical 
color field . The blobs indicate that the corresponding vertices are dressed by the background 
field. 

not exist for quarks, since there is no qqgg elementary vertex. For the same reason, the qqg vertex 
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cannot be dressed by the background field. After some transformations identical to those performed 
in the section 4.3.2, we can represent the NLO contributions to the quark spectrum in the following 
form, which is more appropriate for our purposes: 




(66) 



Compared to eq. (48), two graphs have disappeared, and all the qqg vertices are bare vertices. Like 
in the scalar case, all the NLO terms that involve gluonic vacuum fluctuations (i.e. all the above 
terms except the graphs [f,h]) can be related to the LO spectrum as follows, 



dN„ 



d 3 p 



Jd 3 u [a-T u +(3-T%] 



s 

1 f d?k f 

+ 2^/ (27r)32fc J dZ ^v [a_ fcAa • T u + C-kXa ■ Tg] [a +kXa ■ T v + ( +kXa ■ Tg] 



dN q 
d 3 p 



LO 

(67) 



All the differences between adjoint scalars and quarks are now hidden in the definitions of the 
objects /3 and C±fcAa, that are now given by 



p{u) 




, CfcAa(u) = 



kXd 




(68) 



and in the definition of the operator that replaces T*. This new operator generates shifts of 
the initial condition of a spinor on the Cauchy surface S, and eq. (54) is replaced by 



£-T2 = (&<(«)) 



8{M l (u)) 



(69) 
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